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Abstract 

In the present paper the classical point symmetry analysis is extended from partial differ- 
ential to functional differential equations with functional derivatives. In order to perform the 
group analysis and deal with the functional derivatives we extend the quantities such as in- 
finitesimal transformations, prolongations and invariant solutions. For the sake of example the 
procedure is applied to the continuum limit of the heat equation. The method can further lead 
to significant applications in statistical physics and fluid dynamics. 

1 Introduction 

In the paper we consider such functional differential equations which can be regarded as the extensions 
of partial differential equations. The key idea is that the discrete set of independent variables 
(yi, 2/2) • • • j Un) hi a partial equation is replaced by a continuous set of infinitely many variables 
denoted by [2/(2;)]. In order to illustrate the above extension for partial differential equations to 
functional differential equations we introduce the example 

8f A Of 

where / = f(t; yi, y 2 , . . . , y n ). Taking the continuum limit we obtain 
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where / = f(t; [y(x)]) and the partial derivatives in (0) have been replaced by a functional derivative 
5/5y(x) which can also be denoted by d/dy(x)dx, (cf. Gelfand & Fomin (1963) for the definition of 
the functional derivative). The latter notation is more convenient in some situations (Hopf, 1952). 
The similar procedure of deriving the functional differential equations as the limit case of partial 
differential equations is used also in the work of Hopf (1952) and Breuer & Petruccione (1994). 

Functional equations are used to describe problems of statistical mechanics where the probability 
distributions of phases and their time evolution are studied. In a discrete medium the phase space at a 
given time instant t is defined by velocities and positions of all particles (xi, x 2 , . . . , x n , ui, u 2 , . . . , u n ) 
contained in a considered domain. In a continuum limit (eg. in hydromechanics) the phase becomes a 
continuous function of spatial variables u = u(x). The time evolution of its probability distribution 
is determined by a functional equation. 

Many examples of the functional equations with functional derivatives can be found in physics. 
The areas of interests of the authors are the turbulent flows and turbulent reacting flows. The idea of 
description of turbulence in terms of the characterisctic functional has been introduced in a seminal 
work by E. Hopf (1952). The approach is presented in the book of Monin & Yaglom (1971), we 
mention also the work of Lewis & Kraichnan (1962). The similar functional formulation can also be 
derived for the Burgers equation. Functional formulation of Burgers equation is studied in the works 
of Breuer & Petruccione (1992) and Breuer et al (1996). The authors perform stochastic simulation 
of the Burgers model. The method is based on a discrete master equation which is equivalent to 
the Hopf functional equation in the limit of continuous space. The approach has been continued in 
works of Breuer & Petruccione (1994), Biechele et al. (1999), Friedrich (2002). Another functional 
differential equation is used for the statistical description of a turbulent premixed flame (see Oberlack 
etal. 2001). 

The approach of Hopf originates from the field theory where the functional equations and func- 
tional derivatives are also used (cf. Itzykson & Drouffe 1989, Schweber 1962). The examples of 
functional equations in field theory are the Schwinger equations (cf. Pester et al. 2002) and Wheeler 
de Witt equations (cf. DeWitt 1967, Barvinsky & Kiefer 1998). 

The symmetry analysis based on the Lie group theory has become a powerful tool of analysing, 
simplifying and finding solutions of partial differential equations (cf. Ibragimov 1994, 1995, 1996, 
Cantwell, 2002). The method also gives a deep insight into the underlying physical problems de- 
scribed by the differential equation. Examples of its applications include problems of fluid dynamics, 
where a broad range of invariant solutions for turbulence statistics were found (Oberlack, 1999, 2001). 
However, much less attention have been given so far to the symmetry analysis of functional equa- 
tions. Some of the previous works concern the integro-differential (Zawistowski 2001, Chetverikov & 
Kudryavtsev 1995, Roberts 1985) and delay differential equations (Zawistowski 2002, Tanthanuch & 
Meleshko, 2004). However, to the best of the authors knowledge, the type of functional equations 
considered in the present paper, has not been studied so far in terms of the symmetry analysis. In the 
present work the classical symmetry analysis is extended to functional differential equations which 
contain functional derivatives. The structure of the paper is the following: first, we introduce the 
necessary notation to study functional differential equations. In the main section we extend classical 
Lie group methods to functional differential equations by extending quantities such as infinitesimal 
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transformations, prolongations or invariant solutions. Finally, as an example, the new method is 
applied to a continuum limit of the heat equation. This leads to the transformation groups as well 
as to invariant solutions of the considered functional equation. 



2 Notations 

In the present work we study the extension of the partial differential equation for a scalar function 
. . . , u n , t\, . . . , t m ) of m + n independent variables. The general form of the differential equation 
describing $ writes: 

F( Ul , Un, ti, . . . , t m , $,$,$,...,$) = (3) 

12 q 

where <3> denote the k-th derivatives of the function $ with respect to any possible combination of 

k 

independent variables and q is the highest order of derivative present in Eq. EJ In the considered 
continuum limit u becomes a function of the continuous variable x. For the sake of generality we 
consider that u = u(x) where u = (ui, U2, ■ ■ ■ , u p ) and x = (xi, x%, . . . , x r ) are vectors in a p- and 
r-dimensional space. In this case $ becomes a functional 

$ = $([M 1 (x)],...,[M p (x)],t 1 ,...,t m ) = $([ll(x)],ti,...,* ro ) (4) 

and the partial differential equation (JHJ) becomes a functional differential equation: 

F([ U (x)],t 1 ,...,t m ,$,$,$,...,<l>) = 0; (5) 

12 q 

here again, $ denotes all possible derivatives of order k, which can include partial derivatives with 

k 

respect to tj and functional derivatives with respect to w a (x). The following, equivalent notation will 
be used for the first functional derivatives: 



®,u a ( x) = ~J7~77^ = Q „. TZaJZ - ( 6 ) 



5$ 9$ 
8u a (x) 9w a (x)dx 

The different types of second order derivatives in Eq. will be denoted by 



while the order of any derivative is commutative. Higher order derivatives can be expressed in an 
analogous way. 

As x is a vector in a r-dimensional space, we will use the following, convenient notation for the 
integrals with respect to x: 



G 




clx = / / . . . / dxidx2 . . . dx r (7) 



where G is a domain in r-space to be specified for each particular problem. 

For further purposes we recall here some differentiation and integration rules for functionals. We 
also solve a simple functional equation by the method of characteristics. For the sake of clarity 
we will first present necessary formulae for the partial differential equation © and introduce their 
counterparts in the continuum limit (JSJ. The two approaches will also be called "classical" and 
"continuum formulation" , respectively. 

A function $(tti, . . . , u n , ti, . . . , t m ) of a finite set of variables is differentiable for a particular value 
of its arguments if its variational (or functional) form 5$ is linear in 5ui, i.e. the following relation 
holds (cf. Hopf, 1952) 

» = E^«.+E^.- (") 

i=i i=i 

The analogous formula in the considered continuum limit writes 

A functional derivative of $ exists if its differential form £<3> can be written as the integral (JHJ). 
Let us consider the following function of a finite number of variables $ = YH=i a i u i- Its derivative 
with respect to Uk where 1 < k < n is d^/du^ = J27=i a i$ik — a k where 5ik is a Kronecker delta. We 
consider the following continuum limit of the function $ 



$ = / a(x')u(x')dx' . (10) 
Jg 

The differential form of this functional writes: 

5$ = / a (x')5u(x')dx'. (11) 



G 



Hence, by comparison to Q the functional derivative 5§/5u(x) = a(x) where xeG. Another func- 
tional $ = u(xi) where x% is a given point in the domain G, can also be written in the integral form 

$ = / 5(x — x\)u(x)dx and 5$ = / 5(x — xi)5u(x)dx, (12) 
Jg Jg 

also in this case the functional derivative can be found by comparison with ©, 5&/8u(x) = 5(x — Xi). 
Calculation of second order derivatives is presented on another example below. We consider the 
following function, together with its continuum counterpart 

n n r r 

$ = / j y CjjUjUj — > $ = / / c(x, x')u{x')u(x) dxdx'. (13) 
i=i j=i ' J g Jg 
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Corresponding derivatives of the first and second order write 



9$ ^ V"^ 5$ f / ft /\ / /\ 1 / ft A \ 

^— = > , c fci M i + > J c ikUi -> r p , A = c(x ,x )u(x ) dx (14) 

+ / c(x, x")u(x) dx. 



G 



<9 2 $ <5 2 $ 

Cfcm + C mfc -> t , „ u , , = C(X ,x) + c(x,x'). (15) 



i)iij,()n,,, "' Su(x")5u(x) 

Derivatives of higher orders of the considered function and the corresponding functional are zero. 

We now solve the following hyperbolic equation in a classical and continuum formulation by the 
method of characteristics 

BF >A OF BF f b 5F , , 

9$ dui <9$ 7 a <5w(x') v ' 

where F = F($>, u±, . . . , u n ) in the classical formulation and F = -F($, [w(x)]) in the continuum limit; 
Ui and m(x) constitute sets of independent variables and $ is a dependent variable $ = . . . , «„) 
or $ = $([w(x)]). The characteristic equations of (|16p 

d$ 5$ 

= du\ = ■ ■ • = du n — ► — = 5u(x) for each ie (a, 6) (17) 

determine n integration constants Cj in the classical formulation and an infinite set of integration 
constants C(x) in the continuum formulation. The constants can be employed as new (dependent and 
independent) variables of F. Corresponding solutions of Eqs (fTKJ) have the forms F = F(Ci, . . . , C n ) 
and F = F(Ci, [C(x)]). A few examples of possible solutions of the characteristic system (JT7j) are 
presented below. We can e.g. consider equations 

d$ 

— = dux, du 1 = du 2 , du 2 = du 3 , . . . , dw n __i = du n 



or — = — dui, dui = du 2 , du± = du 3 , . . . , dui = du n 
$ n 

i=l 

to obtain the following integration constants 

$ n * ( ^ 

Gi = or Ci = <Pexp > , 

exp«i \ n *-^t J 

C 2 = u x - u 2 , C 3 = u 2 - u 3 , . . . , C n = n n _i - u n , 
or C 2 =ui- u 2 , C 3 = ui - u 3) . . . , C n = ui- u n . 
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Integration constants (JTSJl have their counterparts in continuum formulation, 



1 



/ 



b 



c 



1 



$ exp [— u(xi)] 



or 



Ci = $ exp 



b — a 



u(x)dx 



C(x) 



du(x) 



dx 



or 



C(x) = — 



(19) 



dx 



where xi is a fixed point in the domain xe (a, b). Hence, a functional that constitutes a solution of 
Eq. (jlfij) in its continuum limit may be written as 



3 Finite and infinitesimal transformations. 

In this section we recall the classical symmetry method which can be used to analyse the partial 
differential equation (jSJ) and present its continuum extension for the functional differential equation 
(J3J). By "symmetry transformation" we understand such transformation of variables which does not 
change the functional form of the considered equation. This means that, for example Eq. (jSJ) in 
the old variables $, t\ . . . t n , u(x) and the same equation written in new, transformed variables 
ti ...i n , u(x) 



are equivalent. Note that x is not transformed in the present approach since it constitutes a contin- 
uous "counting" parameter, such as in a summation for the classical counterpart. Here, we consider 
only such transformations of variables which constitute Lie groups, i.e. they depend on a continuous 
parameter e and satisfy group properties, such as closure, associativity and containing the unitary 
and inverse elements. 

Tabled presents the comparison of finite one-parameter Lie point transformation for the classical 
and continuum formulation. As can be seen, the transformed variables <£>, u(x), t\...i m become 
functionals in the continuum limit and depend on the infinite set of independent variables [u(x)]. It 
should also be noted that instead of the finite set u\ . . . u n , in the continuum formulation we define 
u(x) = (^) x , which is an explicit function of the variable x, since u defines a new variable at each 
point x. This has important consequences in the further considerations. 

For the subsequent purpose of symmetry analysis all variables of equation (J2J), i.e. the sets 
ti, . . . ,t m ,ui, . . . ,u n , as well as $ and all its possible derivatives of any order will be treated as 




(20) 




(21) 
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Table 1: Comparison of one-parameter Lie point transformation for the classical and continuum 

formulation. 

classical formulation continuum formulation 



i. 


$ : 




, U n , fcl) . . 


■ > ; — ) 


$ = [u(x)], 


it. 


Ui 




• ) ^71; ^1 ) • 


• • ; trm 


u(x) = X ($, [u(x')], x, ti, . . . , t m , e) 












xeG 




U n 




• • ) Ufif t\j . 


■ ■ j tffij SJ 




in. 


h ■ 


= tl ($,«i, • • 


■ ■ i ^n; ^1 ) • 


• • ) ^TO) ^) 


ti = tl ($, [u(x)],tx, ...,t m ,e) 








. . . , U ni t\ , 


• • • j ^mi 


= 0* m ($, [u(x)],ti, ...,t m ,e) 



independent variables. Now, the following, new differential operators are introduced: 



V d ^ d d d . . 



j=i - j j=i 



= 1, . . . , m 



j=i • j j=i 



i = 1, . . . , n. 



The partial derivatives e.g. of the form d/dti will act only on terms which depends explicitly on tj. 
Within this formulation the derivatives of $ can be expressed as: 



Analogous definitions will apply in the continuum limit (jHJ). The derivatives (J22j) have the following 
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counterparts 



V d d d 

— = 1 '-'-7^ ■ — 



Vt k dt k ,fe <9$ ^ >^ 3 d§ t . 



+ £ / dx$ 



a=l 



,tfcUa(x) 



I? 5 ^ d d . . 

+ + 2^^^)t J T^~ (26) 



p r 5 



(for a = 1, xeG) and the derivatives of $ in terms of the new differential operators write: 



< ^ > ,U a (x)l/ /3 (x / ) 



£>t/ £>u a (x)dx' 



2>u Q (x) dx Vu p (x' ) dx' ' 

'" a(x)ij = Dii Q (x)dx^' ^ tk = WjVh' ( 2) 

Note that all derivatives in (J2?j) are commutative. Also it is important to distinguish between x and 
x' which denote different integration indices such as i and j in two consecutive summations. The 
quantities given by formulae (i), (ii), (Hi) and derivatives of $ can be written in a Taylor series 
expansion about e = 0. Their infinitesimal forms, after neglecting terms of order 0(e 2 ) are given in 
Table El A few notation particularities should be noted. Indices of ( are separated by a semicolon 
to distinguish it from derivatives. Functional £ x which denote infinitesimals corresponding to u is 
an explicit function of x. The same dependence holds true for infinitesimals corresponding to the 
functional derivatives of $, such as C,u(x)- in these cases the index of the set [u(x)] has been given a 
different name such as [u(x')] to avoid confusion with the parameter x. The remaining infinitesimals 
do not depend on x. 

The key property of Lie group method is that the finite transformations, given by the formulae 
(J,)— (Hi) can be computed from their infinitesimal forms (iv)-(vi) (cf. Bluman and Kumei, 1989). 
According to the Lie's first theorem, the finite form of the transformation can be obtained by in- 
tegrating the first order system of equations. For the continuum formulation this system takes the 
form: 

d$ dti du a Cx) . . 

— = 77, — = % = 1, . . .,m, dg = U, a = 1, . . . ,p, xeG; (28) 

where the latter equations should be integrated with the initial condition 

e = 0: <5 = $, ii = U, u a (x)=«a(x). (29) 
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Table 2: Comparison of infinitesimal transformations for the classical and continuum formulation. 



classical formulation continuum formulation 



iv. 
v. 


$ = 

Ml = 


$ + r/($,Mi, ...,M n ,ti, ...,t m ) £ 

= Ul +^i($,Ui,...,14 n ,*i,...,t m ) e 


$ = $ + 77($,[u(x)],fi,...,t m )e 
u(x) = u(x) 




vi. 


U n = 
tl = 


: u n + £„(<!>, tii, ...,u n ,ii,...,i m ) £ 
*i +6i(^,Mi,...,M ra ,ti,...,t m )e 


+£ x ($,[u(x')],x,i 1 ,...,i m ) 

ti = ti + f tl ($> [u(x)],t!,...,t m ) 


e 
e 


vii. 


tin 


: + & m ($,til, ...,M„,tl, ...,t m )e 
= $,t! +C;ti(^>Wl>->Wn>*l>->*m)e 


*m=tm + ^t m (^,[u(x)],*i,...,i m )e 

= $ tl +0 1 ($,[u(x')],ti,...,t m )e 


viii. 


d> _ 


= ^.t™ +C;tm(^,«l,-,«n,*l,-,*m)e 
^,1*1 "I - C;ui (^*, til, tin, ^1, ^m) £ 


?,t m =^,t m +C;t m (^[u(x')],* 1 , 


tm) £ 






= ®,u n + C,u n (®, til, -, ti„, *1, -, <m) £ 


+C;u a (x)(^, [u(x')], X, t 1} t m ) £ 




ix. 




$,« Q (x)« /3 (x') = $,« a (x)« /3 (x') 








+ C;u^($, w l, — ,«n,*l) — ,*m)£ 


+C;MxW(x') ( $ > [ u ( x ")] , x , x ', tl , 


•••■> ^m) ^ 
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Now, our aim is to find the infinitesimals forms rj, £ QX . To do this we should first express the 
infinitesimals ( in terms of i], £ t ., £q, x and independent variables t±, . . . , t m , [u(x)], $. In the classical 
formulation the following relation is used for this purpose 



T>Ui 



k=l 



k=l 



n 



fc=i 



in 



T>4>, 



v^m~ k ^ Vui vt k ~^^' Uk Vui T ^ T,tv p« 



(30) 



fc=l 



its continuum counterpart writes 



X>w Q (x)dx 



(3=1 



E/4 



,u,a(x')" 



^ x ,($,[u(x)],x',t,£) dx/ 



£>u a (x)dx ' fe £>u Q (x)dx 

When the infinitesimal forms (iv)-(viii) , are introduced into equation (J3TJ) we obtain: 



(31) 



l>u a (x)dx 



E / ( $ ,«M X ') + C;«M x ') e 
0=1 ^ G 



D(Mx')+^) dx , 



2>u Q (x)dx 



+ E($,t fc + C ; t fc e 



fc=i 



V(t k + £ tk e) 
Vu a (x)dx. 



(32) 



Equation (|32|) can be further split into two equations, containing terms 0(1) and 0(e), respectively. 
The first of the two gives the identity 



2>u a (x)dx 



Ef ^ Vur(x!) , , 



/3=1 
P 



"D« a (x)dx 



E <W / $,^( X ')5(X - x') dx = $ Uq 
/3=1 ^ G 



(x) 



(33) 



where Dup{x!) /Vu a (it)&x. = 5 a p5(x. — x') has been used where 5 a p and 5(x — x') = 5(xi — x' x ) • 
• ■ • 5(x r —x' r ) denote the Kronecker and Dirac delta respectively. From 0(e) we obtain a formula for 
the infinitesimals C;u a (x) 



Vr] 



2>u Q (x)dx 



= E 

(3=1 
in 



^^x' . > Vup(x! 

+ W(x'); 



Dw a (x)dx 



2>u a (x)dx 



dx' 



fe=l 



2>u a (x)dx' 



hence, 



C;«a(x) 



Vr] 
Vu a {x)dx 



0=1 



E / *.«, 



(x') 



2}u a (x)dx 



fc=l 



£>w Q (x)dx 



(34) 



(35) 
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By analogy, formula for the infinitesimals (. tl can be found 

C*- = ^-E I * M*0^dx'-f>^. (36) 



0=1 fe=l 
The infinitesimals of higher orders will follow from the following, recursive formulae: 

^,n a (x(i)),...,n )3 (x(-i)) _ ^ f A ^ 

5=1 



; ,n a (xW),...,n /3 (x("-D) \^ [ & ^x , 

D % (x«)dx« * (xW) ^ (xM)) ' U5(x) % 7 (x^))dxW 

£> % (x«)dx« 



+ $ ,« a (x(i)),...,« (3 (x(-i)),t fc — • ( 37 ) 



or 



^,« Q (x(i)),...,n /3 (x(»-D) y * ^ 



2^ / *,«a(xW),...,U P (x(-l)),« 4 (x)-^— OX 
js_i ^ G 



+ X^^« a (xW),..., U/ 3(x(-i)) A ^ i - (38) 
fe=i •? 

4 Generator X and its prolongations. 

Once all the necessary infinitesimal forms are obtained they can be substituted into the equations 
@ or (J3J) written in the transformed variables. In order to simplify notation we will assume that 
Eqs (JHJ) and (jSJ) only contain derivatives up to the second order. The generalization of the following 
relations to the case of higher order derivations is straightforward. After expansion in Taylor series 
about e = in both, classical and continuum formulation the expanded equation has the form: 

F + eX^F+j [X^] 2 F + 0(e 3 ) = (39) 

where X^ 2 ^ in the classical formulation is given by the formula 

d m d n d m d n d 

x(2) = ^ + E^flT + + 



J l J J=l J 3=1 ' J 3=1 ' 3 



mm q n m q n n q 

q^— + Cw*» q^— + w ( 4 °) 
j=i fe=i j=i ^=1 ^ j=i fc=i ^ 



and is called the prolongation of the generator X 



3=1 J 3=1 J 
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of the second order. The corresponding formulae for the continuum limit write 



f) m f) p r 

7=1 J a=l J( ^ 



5 

ax 



(42) 



and 



X (2) 



P r mm q 

+ E / dx u + E E fc* ^— 

p m „ 

+ EE / dx c^« 

a=l fc=l >/G 



uv ,« Q (x)i t 



p p r r § 
+ EE / / dxdx/ C;« Q (x)^(x')T^ • (43) 



a=l /3=1 



'GiG ^,« c «(x)w / g(x') 



The first term in Eq. (|39j) equals zero, as follows from the equations (|3*j) or (J5jl . All the remaining 
terms [X( 2 )] n in ()39j) . representing a successive application of X^ 2 ) will be zero if the following relation 
holds 

X (2) F = 0. (44) 
In order to find the infinitesimal transformations we use the condition 

[X( 2 )F]| F=0 = (45) 

where, in the continuum formulation, the prolongation X 1 - 2 ) is expressed by the formula (|4*3*|) and the 
forms of infinitesimals ( are found from relations (|33|) - (j38)) . The resulting condition constitutes an 
overdetermined system of linear differential equations. In the continuum limit we obtain a set of 
functional differential equations. This system can be further solved for the infinitesimals T], £ f . and 

£ax- 

5 Invariant solutions 

If the functional differential equation (J5J) admits a symmetry given by the generator (|42jh then a 
solution $ = 0(ii, . . . ,t m , [u(x)}) of this equation is called an invariant solution if it satisfies the 
relation 

X[$-e(t 1 ,...,t m ,[u(x)))) = 0. (46) 
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After employing ()42j) and expanding the derivatives, from (jlojl the following, hyperbolic functional 
equation is obtained 

i=l a=l v ' 

This equation can be solved by the method of characteristics. The corresponding system of equations 
writes 

— — = • • ■ = — — = — = — ^ for each x e G and a = 1, . . . ,p. (48) 

6l 6m V Sax 

Above, has been replaced by Note that the last term in fact corresponds to an infinite set of 
equations for each a and each point in G. The infinite set of constants, which is a solution of the 
above system, can be employed as new variables in Eq. (|5f . As in the considered case one of them, 
say Ci, will be a dependent variable and the rest will constitute a set of independent variables, the 
following relation holds 

C x = H(C 2 , . . . , C m , [d(x)], . . . , [C p (x)]). (49) 

After the process of solving characteristic system for a partial differential equation with a finite set of 
variables, the number of independent variables is reduced by one. In the case of functional differential 
equations, in formula (|4*H|) one point of the considered domain will be excluded from a set [Cj(x)], 
xe G/ixjt}. 



6 Example 

For the sake of clarity we will consider, as an example, the continuum limit of a heat equation in 
infinite many dimension. In a classical formulation we take m = 1. Hence, the equation for a function 
of n + 1 variables $ = $(«!,..., u n , t) writes 

^ = ^^?' (50) 
i=i 



while in the continuum limit we consider 



m = J a du(x)d x du( x )dx dx ° r ^ = I ( 51 ) 
where $ = <&([u(x)],t). The generator X of Eq. (|HT|) (cf. Eq. (j41jl ) has the form 

d d f b 5 
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The second prolongation X^ 2 ^ of X necessary for (JoTj) (cf. Eq. (jjOj) ) is defined as 

d 



x,2> - "4 + 4 + / di ^ +C;< 



9$ 



+ / dx Ciu(aOii(*) , (53) 

Ja ,u(x)u(x) 

where any unneeded ( has been omitted. Applying (J53J) to (JoTj) the first three terms of (Jo"Tj) have no 
effect, while the fourth one acts on $ t to lead to (. t . The last term in (p)3*|) . acting on (JoTj) can be 
written as 

b ^ rb 

dx(; U ( x)u ( x) — / ® jU ( x >) u c x >)dx' = (54) 



dx(.,u(x)u(x) / 5(x-x')dx = / C^^mW dx. 
As a result we obtain 

C;t - / C^n*) dx = 0. (55) 

Into the above equation we substitute the infinitesimals (-t, C;u(x)u(x), found from Eqs (fHKj) - (fHTj) . Their 
forms, without derivation, are given below 



T>u(x)dx T>u(x)dx T>u{x)dx J a T>u{x)dx 



b 



^.utrHi') ^ / \ i dx — $, u (x)t; 



T , M(x)dx ' Vu(x)dx 



Vu(x)dx \ ' Vu(x)dx J 
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If the definitions of differential operators (J2fijl are used in the above equations one obtains 

^ dt •* at >*a$ 

- jf ( $ .«W^ + *,t.(-)*,*^) d * and 

£«(*)«(*) = ^^ + 24) ^)^y^ + Cl) i n W ^2 «(*)«(*) 

" A $ '" ( *' } + 2$ ' uW feR¥ + *.«<«) W + ^ijjj d * 

'* V«J«(x) 2 ,tt(i8) <fa(z) 9$ + + ^"^^ j 

" A + arj d * - ^ + $ m.) M J ■ 

Next, we use Eq. (J5T|) in order to substitute for the derivative $ t 

$ t = / ds - 

a 

The resulting equation has the form 

d ^-/ -5T + 2 ^,,^^ - / A^2 dx dx 



<9£ J a 5u(x) 2 J a ,u ^ x ' \ dt d&du(x)dx J a 5u(x 

rb s2c ac \ rb rb 



+ I (/ ^^-f ) dx + 2 / / $ ^o^ d ^ 

2 L L ^^^ dxdx ' + 2 J a I W)^) ^ )dj d * d ^ 



/ / $2 n(x) $ ,«M^| :da;d:E, + / / ®M*>(^®Mx)-Q^I dxdx ' 

*J (I %f CI. tJ Ob *J CL 

+ 2 1 ^^k dx+2 L ^)^)§ d -=o. 



15 



As the terms rj, £ t , £ x do not depend on the derivatives of from Eq. (|58|l the following system of 

differential equations can be obtained, where on the left hand side the coefficient function is written: 

: H = °> ( 59 ) 

: = 0, for each are (a, 6), (60) 

^•Wut^^W : = 0, for each xe(a,b), (61) 

* 2 .W ■■ = 0. (62) 

«U*M : I' ^ d*- - § + 2^ = 0, (63) 

for each x e (a, b), 



for each x, x' e (a, 6) , and x ^ x', 

,«(*) • + z «9$<9w(x) dx J - 5«(x') 2 ~ ' 1 j 

for each x e (a, 6) 

Both equations (|6l3|) and follow from the second line of Eq. ()58jl . To see that this is true we will 
differentiate Eq. (f5SJ) with respect to by 5/<^,«(a:")u(*"')' Since $ is an explicit function of two 
variables it follows that 

6(x - x", x' - x'") = 8(x - x") 5{x' - x'"). (67) 



,u(x)u(x') 



S^,u(x")u(x"') 

On the other hand, $> tU (x)u(x) is a function of variable x only, say $ jU (x)u(x) = 9( x )i hence 

§^^ = ^L = 5( X -^). (68) 
o® ,u(x")u(x") og[x ) 

The functional derivative of $, u r x )u(x) with respect to & jU {x")u{x'") is zero if x" ^ x'" . Finally, after 
differentiation of (J55J) with respect to $, u (x")u{x'"), f° r x " — x "' we obtain 

r^-^dx^^dx'-f) (69) 

+2 f a 6 / a 6 5(x - x")8(x' - x")^krdx dx' = 
whereas in case when x" ^ x'" we have 

2 f [ 8(x-x")6(x' - x'")4^r + 5{x' -x")5(x-x"')J$-dxdx' = 0. 
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Note that the presence of two terms inside the integral is a consequence of an equality & >u (x)u(x') — 
$ )U (x')u(x) ■ From the above equations formulae and (JMj) are clearly obtained. 

The above system of equations is now solved in order to find the form of infinitesimals rj, £ t and 
£ x . From relations (|59j) - (|fi2jl we obtain: 

Ct = (70) 

C = UW*%x,t), (71) 

V = /i(Kz%*)$ + / 2 (Kz')],f), (72) 

Next, we consider formula (j63j). We recall that £ x is an explicit function of x, hence, the general 
solution of (jHSJ) has the form 

£■ = ^1 [ u(x')dx' + H(x,t,[u(x')))+f 3 (x,t) (73) 

where H(x,t, [u(x')]) is a functional depending explicitly on x, such that its functional derivative 
6H(x,t, [u(x')]) I Su(x) = 0. Without an explicit dependence of £ x on x the second term on the RHS 
would disappear as we would obtain H(t, [u(x')]) = 0. Now, the form (J73|) is introduced into Eq. (JHlj) 
to give 

lrffr , SH(x,t,[u(x")}) ldCt 5H(x',t,[u(x")]) , , ( . 

2 dt + ^(x') ~ 2 dt ' 1 J 

Below, we prove that the functional derivative of H(x, t, [u(x')]) with respect to u(x') is a function of 
three variables G(x,x',t). First, we differentiate Eq. (J73j) with respect to u{x') and use the property 
5H(x',t, [u{x")])/5u{x') = to obtain 

Let us now differentiate Eq. (|74|) with respect to u(x\), such that xi / i / x'. We obtain 

S 2 H(x, t, [u{x")\) 5 2 H(x', t, [u{x")\) 



Su(x') 6u(xi) 5u(x) 5u{x\ 



(76) 



Now, if we note that Eq. (J73J) can be written also for pairs of variables Xi, x and x±, x', the LHS of 
the above formula equals 



5 {5H{x,t,[u{x")})\ 5 f5H{x u t,[u{x")}) 



Su(x r ) \ Su(xi) ) 8u{x') \ 5u{x) 
and the RHS is 

5 {5H(x',t,[u(x")])\ 5 fSH{x u t, [u(x")]) 



5u(x) \ Su(xi) J 5u(x) \ 5u(x') 
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(77) 



(n 



which is an expression with opposite sign to that in ()77|) . It follows that the second order functional 
derivative of H is always zero, hence 

From Eq. (|75j) we obtain the relation G(x,x',t) = —G(x',x,t) — d£ t /dt, additionaly G(x,x,t) = 0, 
as the first derivative of H(x,t, [u(x")]) with respect to u(x) is zero. Now, the formula (JTHj) for the 
infinitesimal £ x writes 



^ x = f u ( x ') ^ x> + f G(x,x',t)u(x')dx' + f 3 (x,t) 

The above form can be rearranged as follows 

C(x, a/, t)u(x) dx' + /s(x, i) 



*0) 



where 

C(x,x',t) = -C(x',x,t). for x^x' (82) 
C(x,x,t) = (83) 
This result, together with (J72j) can be further substituted into to give 
(J/i(K^)],t) 1 /- fe aC(a;,x',t) . /NJ , l9/ 3 



, ufzW ^. (84) 

5u(x) 2 J a dt K ' 2 dt K J 

After integration we obtain 

Mu(x%t) = -Uala dC %f' t] u{x')u{x)dx'dx (85) 

~\ f^u(x')dx' + U(t). 

If we recall the differentiation rules presented in Section El we find that the functional derivative of 
the first RHS term in the above equation is a sum of two integrals (cf. Eq. ()14Jl ) and is equal to 
the corresponding term in Eq. (|8%)l only if for x 7^ x', dC(x,x' ,t)/dt = dC(x', x, t)/dt. It further 
follows from (JH2*j) that the only possibility is dC(x, x', t)/dt = for x 7^ x'. Finally, the form (|72*|) is 
substituted into equation (jfifij) to give 

$ (-l d ^I b a l"c(x,x'Hx)u(x')dxdx' 
-|L 6 ^^)d^ + f + l^(6-a)) 
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where c(x, x) — 1 and c(x, x') = for x 7^ x'. As £t, /3 and f± neither depend on $ nor on u(x), from 
the above equation it follows that 



d% 
dt 3 

d 2 h(x } t) 
dt 2 



(87) 
(88) 



dfi (b — a) d 2 £ f , , 

df 2 f b S 2 f2 



, . 9 dx = 0. (90) 
dt J a 5u(x) 2 v ; 

Hence, the final solution of the system (J59j) - (j66j) writes 

C t = a l t 2 + a 2 t + az (91) 

= / C(x,x',t)u(x') dx' + a 4 (x)t + a 5 (x) (92) 

J a 

i] = — \~^ a i J J c(x, x')u(x)u(x') dx dx'+ 

dx+ ^ - ^- ait + aA $ + f 2 ([u(x)],t) (93) 

where C(x, x,t) = |(2ait + a2), C(x,x',t) = C(x,x f ) = —C(x',x), and fi is a solution of equation 

The infinitesimals (l9"T|) - (j9"3)) may directly be compared to those of the classical heat equation (JHUj) 
(cf. Ibragimov 1995) 

& = oit 2 + a 2 t + o 3 (94) 
1 n 

& = -(2a 1 t + a 2 ) Ui + ^2 QjUj + a i4 t + a i5 , i=l,. . . ,n (95) 

3=1 



f 1 - 1 " 1 \ 

7?= - I ^ a i^ M » 2 + o^Qi4Mi + ^nait + a 6 j $ + / 2 (Mi,...,M n ,t) (96) 

\ i=2 i=l ~ J 

where Cy = — Cji for z 7^ j and Ca = 0. As it is seen, the form of infinitesimal £t, Eq. (|94j) . is the 
same as in the continuum limit. In formula (|9"K|) the two first terms on the RHS can be rearranged 
as l/2(2ait + 02) Yli=i C'ijUj with C'^ = —C'^ for i 7^ j and C' u = 1. Then, in the continuum limit we 
obtain the first RHS term of ([92)1 . The constants a i4 and become in (|9"2"j) functions of a continuous 
parameter 2. The double integral in Eq. ()93|) is a continuum counterpart of the first RHS term in 
). We note that the second order derivatives of these terms, present in Eq. (j66|) and its classical 
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counterpart, are the same in both formulations 



dfx 
duk 


1 


1 


d 2 h 
du\ 


1 




Sfi 
Su(x) 


i r b 


(x, x') dx' 


5 2 /i 
5u(x) 2 


1 

Oi. 

2 





The second RHS term in Eq. (|96|) becomes an integral in Eq. (|93J). In the third RHS term, instead of 
the number of variables n, in the continuum formulation (jfJHJ) we obtain the length of the considered 
integration domain which equals (b — a). 

From Eq. (|91j ) -([93 )) we can distinguish the following symmetry groups, admitted by the considered 
Eq. flSIJ) 

X, = |, (97) 

X u ( x ) = g -, r , for each xe(a,b), (98) 

O'u(x) 

^ (5 

X u( .h, = „ (a 0^-«(s)^, (99) 

for each x, x' e (a, 6) and x ^ x', 
5 d 

X u (a;)$ = 2t — u(x)&——, for each xe(a, b), (100) 



d 5 

X tu ( x ) = 2t— + u ( x ) $ u ( x y for each xe(a,b), (101) 



X u ( x ) t $ = t 2 — + fr»(x) - - / / c(x,x')u(x)m(x') dxdx' 

' ? ' ~ ^"^tIL for each x e ( a ' 6 )- ( 103 ) 



2 <9$ 

We will investigate the invariant solutions of Eq. (}5T| under the combined symmetries (J98|) . ()100|) 
for each xe(a, 5) and (jl()2|) . The invariant solution of Eq. (|51|) satisfies the following, hyperbolic 
equation (cf. formula (jUj) ) 

Z" 6 5$ 1 f b 

/ (a 4 (x)t + a 5 (x)) dx = -a 6 $ - -$ / a 4 (x )w(x )dr (104) 

Ja 0M(Xj 2 i / a 
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the characteristic system writes 
dt 
" 



Su(x) 



-a 6 $ - ±$ £ a 4 (x')u(x')dx' <H(x)t + a 5 (x) 
We will solve the following equations 

dt = 0, 

Su(x) 5u(Xl) rr;, 

— 7 — r = — 7 : ; r, for hxed X\ 

a<i{x)t + a<s{x) a 4 (xi)t + a 5 (xi) 

d$ ll Su ( x ) ( 2fl 6 + / tt 6 a 4 (x')w(a; / )dx'J dx 



2 J" (a 4 (x)t + 05(2;)) dx 



which give a set of integration constants 
C(x) 



U[X) 



u(xi) 



t = r, 

for each x e (a, 6), 



a 4 (x)t + a 5 (x) a 4 (xi)t + a 5 (xi) 

4ok f 6 'u(x)dx + f fe f b a 4 (x'WxWx')dx'dx 
$ = C$exp Ja v ; JaJa 4 V ; v ; v ; 



(105) 

(106) 
(107) 

(108) 

(109) 
(110) 

(111) 



4 f (a 4 (x)r + a 5 (x)) dx 

The constants constitute a new set of variables. We treat C$ as a new dependent variable, hence, 
from a relation C$ = /i(r, [C(x)]) we find a formula for $ 

4a 6 J a b u(x)dx + f Q b J a b a 4i (x')u(x)u(x')dx'dx 



$ = ^(r, [C(ar)])exp 



4 J (a 4 (x)r + a 5 (x)) dx 

= h(T,[C(z)])exp{H) 



(112) 



where for brevity the term in the exponent is denoted by H. This result can be further substituted 
into Eq. ()51|) to find the final form of the invariant solution. First, we compute the necessary 
derivatives in terms of new variables 



<9$ dh , N \' b aAx)dx 

— = — exp (H) + hexp (H)H— h Ja V ; , 

dt dr yy ' yy ' J b a ( ai (x)T + a 5 (x))dx 



(113) 



the expression for the second order functional derivative of $ is lengthy, however, after integration 
with respect to x all terms that contain functional derivatives of h will cancel. Below, we write the 
final form, after integration 

rb 1 2 



rdx 



h exp (H) 
h exp (if) 



2a 6 + f a4(x')u(x')dx' 
2 J a 6 (a 4 (x)r + a 5 (x)) dx 
f fl 6 a 4 (x)dx 
2 J a b (a 4 (x)r + 0.5(0;)) dx 



(114) 
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We can now compare (jll3|) with (jll4|) and perform the time integration to obtain the form of 
functional h. Finally, the solution of Eq. (pTTj) . invariant under the considered symmetries writes 



*(t, Hx)]) 



(115) 



g([c(x)]) 



exp 




V 1 fa 04 ( X ) dx +fa fl 5 ( X ) d: 



IX 



4 J a ai(x)dx(t f a a 4 (x)dx+ a 5 (x)dx^ 



7 Conclusions and perspectives. 

In the present work the classical, point symmetry group analysis is extended from partial differential 
equations to their counterparts in the continuum limit. In particular, we introduce the procedure of 
applying symmetry analysis to the case when functional derivatives are present in the equation. As 
example the method is further applied to the continuum limit of a heat equation and the Lie point 
symmetries, admitted by this equations are derived. From the infinitesimal transformations one can 
also find the invariant solutions of the considered equations. 

The presented extension of the Lie groups can be a useful tool for analysing the functional 
equations. Though we have only given the heat equation as an application of the method we believe 
that the new approach is highly relevant to a variety of important functional differential equations 
(FDE) in mathematical physics. Generally speaking, very little is known on how to analytically treat 
and solve FDEs (numerical treatment is difficult anyway because of the high dimensionality). Hence, 
the methods may give a chance to treat equations which so far have been put aside because of the 
missing analytical methods. In fact, the benefit is twofold since the symmetries not only allow for 
analytical solutions but are also useful in itself since symmetries illuminate the axiomatic properties 
of the physical model equations. 
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